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I.  RESIDUAL  STRESS  ASSESSMENT 


The  existence  of  the  residual  stresses  in  thin-film  coated  structures  has 
attracted  the  interest  of  many  investigators.  The  study  of  the  residual 
stress  behavior  on  evaporated  films  was  first  performed  by  Stoney^  as  early 
as  1909.  He  derived  a  simplified,  although  inaccurate,  expression  for  the 
stress  in  a  thin  film  deposited  on  a  substrate  strip  based  on  the  elastic 
theory  of  simple  beams.  Since  then,  more  work  has  been  done  by  Hoffman,^ 
Finegan,^)  Brenner  and  Sanderoff^  to  refine  the  stress  prediction  tech¬ 
niques  . 

However,  the  residual  stress  prediction  technology  depends  on  accurate 
measurements  of  the  geometric  deformation  of  the  thin-film  coated  struc¬ 
tures.  The  film  stress  can  be  tensile  or  compressive  depending  on  the 
mechanical  properties  of  the  films  and  substrates  such  as  the  Young's  moduli, 
magnitude  of  the  mismatch  of  the  thermal  expansion  between  film  and  substrate, 
the  processing  parameters  such  as  deposition  rate,  film  thickness,  substrate 
temperature  and  chamber  pressure  during  deposition,  and  the  coating/substrate 
adhesion  characteristics.  Deformation  of  the  structures  such  as  change  of 
local  radius  of  curvature,  strains,  deflections,  etc.  will  be  induced  due  to 
the  residual  stresses  existing  in  the  structure. 

Tables  1  and  2,  taken  from  Ref.  2,  summarize  the  currently  available 
stress  data  for  dielectric  and  metallic  films  on  a  variety  of  substrates. 
These  data  provide  some  insight  of  both  the  residual  stress  behavior  and  their 
magnitudes  which  may  be  expected  in  a  thin-film  coating.  Furthermore,  the 
mechanical  properties  of  the  thin-film  coating  dictate  the  type  (tension  or 
compression)  and  the  magnitude  of  the  residual  stresses  in  a  thin  film. 
Coatings  often  exhibit  properties  which  may  differ  significantly  from  their 
bulk  properties.  Hence,  one  would  need  to  measure  the  mechanical  properties 
of  coating  materials  in  the  thin-film  configuration.  The  properties  needed 
include  Young's  modulus,  shear  modulus,  free  thermal  expansion,  and  thermal 
conductivity  as  functions  of  temperature. 
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Table  1.  Intrinsic  Stresses  in  Dielectric  Films  Approximately  5000  k  Thick 

(from  Ref.  2) 


;>y  - 

K 

Substrate  Temp. 

Substrate 

Stress 

y£  Material 

(#c)« 

Material 

(lO^dyne/cm^)^ 

Signc  Method0 

A,  thermally  floating  at  ambient  temperature. 

"Values  in  parentheses  are  relative, 
and  T,  compression  or  tension. 

^B,  end-supported  beam;  C,  cantilever  beam;  and  ED,  electron-diffraction 
technique.  , 

eAl  and  Mg  thalocyanine.  0 


Table  2.  Intrinsic  Stresses  in  Metal  Films  Approximately  1000  A  Thick 

(From  Ref.  2) 


Glass 

0.75 

T 

C 

Copper 

1.0 

T 

T 

Mica 

0.2 

T 

C 

Copper 

0.75 

T 

C 

Cellulose 

1.2 

T 

ED 

Copper 

0.1 

T 

C 

Cellulose 

4.6 

T 

ED 

Quartz 

2.9 

T 

P 

Copper 

0.85 

T 

C 

Copper 

0.9 

T 

C 

Mica 

0.2 

T 

B 

Cellulose 

4.4 

C 

ED 

Copper 

1.5 

T 

C 

Permalloy 


Glass , Silica 


Glass .Mica 


Nickel 


Silica 


A  means  thermally  floating  at  ambient  temperature, 
refers  to  tension,  C  to  compression. 

B,  beam  supported  on  both  ends;  C,  cantilever  beam;  ED,  electron  diffraction; 
and  FR,  ferromagnetic  resonance. 

“Poisson-corrected. 

ePoi88on-corrected,  36.5*  angle  of  incidence. 
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Failure  or  degradation  of  thin  films  as  a  result  of  either  the  deposition 
process,  or  under  external  mechanical  and  thermal  loadings,  requires  a 
detailed  and  systematic  investigation.  Failures  may  include  local  cracking, 
wrinkling,  spallation,  melting,  and  sublimation.  These  failures  can  be  gener¬ 
ated  by  thermal  stress  or  mechanical  stress,  or  a  combination  of  the  two. 
Detailed  investigations  are  hence  needed  to  identify  the  failure  modes  and 
failure  mechanisms  in  order  to  provide  guidelines  to  enhance  the  predictabil¬ 
ity  of  the  response  of  thin  coatings  under  given  external  mechanical,  laser, 
or  nuclear  environments. 

Therefore,  the  activities  of  the  residual  stress  assessment  task  of  the 
program  include: 

1.  Analytical  derivation  of  stress  behavior  in  thin  film  coated  sub¬ 
strates 

2.  Development  of  viable  measuring  techniques  to  determine  the  deforma¬ 
tion  of  thin  film  coated  structures 

3.  Establishment  of  a  testing  methodology  to  quantify  coating  material 
properties  in  thin  film  configurations 

4.  Identification  of  failure  modes  and  failure  mechanisms  to  provide 
information  toward  the  establishment  of  failure  criteria  of  thin 
film  coated  structures 

The  results  and  the  progress  of  each  of  the  above  activities  are  dis¬ 
cussed  in  the  following  sections. 


A.  RESIDUAL  STRESS  IN  THIN  FILMS 

Stoney's  original  problem  was  re-evaluated  to  clarify  the  underlying 
assumptions  and  to  identify  problem  areas  that  were  needed  for  thorough  under¬ 
standing  of  residual  stress  distribution  in  films  of  thickness  t.  The  model 
includes  a  planar  substrate  of  thickness  d  and  a  typical  linear 
dimension  4.  During  the  course  of  the  film  deposition,  the  substrate 
undergoes  deformation.  The  deformation  is  composed  of  two  sources:  thermally 


induced  and  non-thermally  induced.  The  former  is  caused  by  the  temperature 
rise  and  the  temperature  gradient  associated  with  the  film  deposition,  while 
the  latter  is  caused  by  the  atomic  interaction  between  the  film  and  substrate 
materials.  At  post  vacuum,  the  temperature  of  the  film-substrate  assembly 
returns  to  ambient  conditions  and  the  thermally  induced  deformation 
disappears.  The  deformation  of  the  system  is  therefore  due  to  the  "non- 
thermally  induced"  stress  as  a  result  of  film  deposition.  This  stress  is 
defined  as  "intrinsic  stress".  Therefore,  the  intrinsic  stress  is  regarded  as 
the  level  of  stress  which  exists  in  the  system  as  a  result  of  the  film 
deposition  without  any  boundary  constraint  to  the  system.  Therefore,  any 
thermal  gradient-induced  stresses  and  stresses  associated  with  boundary 
condition  changes  are  not  included  in  the  intrinsic  stress.  The  general  case 
as  well  as  those  for  both  cantilever  strips  and  circular  disks  are  considered. 

The  following  assumptions  are  adopted: 

1.  The  substrate  material  is  homogeneous,  isotropic  and  linearly 
elastic  with  Young's  modulus  E  and  Poisson's  ratio  v. 

2.  t  d  <<  l 

3.  The  edge  of  the  assembly  is  traction-free. 

4.  The  bond  between  the  film  and  the  substrate  is  perfect. 


5.  The  film  experiences  an  average  two-dimensional  hydrostatic  state  of 
stress  which  is  uniform  throughout  the  film  thickness.  This  average 
instrinsic  stress  is  denoted  by  S. 


/  -r 
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Take  cartesian  coordinates  x,  y,  z  with  the  plane  z  -  0  coinciding  with 

the  central  plane  of  the  substrate.  Therefore,  the  coordinate  2  runs  across 

the  thickness  direction  of  the  substrate.  These  assumptions  apply  to  the  case 

where  the  substrate  is  subjected  to  a  state  of  spherical  bending.  Thus,  one 

has  the  in-plane  linear  stress  components  o  (z),  0  (z)  as 

xx  yy 

°xx(z)  =  °yy(z)  5  °{z)  ( 1} 

=  0  otherwise  (i,j:x,y,z) 

and  the  in-plane  linear  strains  e  (z)  and  e  (z)  obey 

xx  yy 

£xx(z)  =  eyy(z)  =  e(z)  =  e(z)  =  '  2~rr^  z  o(z)  (2) 

where  r  is  the  radius  of  curvature  of  the  neutral  plane  and  c  is  the  distance 
of  the  neutral  plane  from  the  plane  z  =  0.  From  Eq.  (2)  the  stress  in  the 
substrate  can  be  expressed  as 


o(z) 


E  c  t  z 
1  -  v  r 


The  requirement  of  force  and  moment  equilibrium  over  the  cross  section 
prov ides 


k'- 


> . 


Is 


d/2  +t  d/2  r  „ 

J  Sdz  +  J  o{z)dz  =  St  -  /.■  — 7-  =  0 
d/2  -d/2  U  v;r 


d/2  +t  d/2 


cd3E 


J  Szdz  +  J  o(z)zdz  =  x  t(t  +  d)  -  .  — r—  =  0 

d/2  -d/2  d  U  v;r 


The  film  stress  S  and  the  distance  c  between  the  neutral  and  middle  plane  are 
therefore  solved  as 
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Denote  the  displacement  in  z-direction  of  a  point  on  the  neutral  plane 
by  u>  and  let 


3w  Bui  _  q 
3x  "  3y 


at  x  =  y  =  0. 


Since  the  deformation  w  is  axisymmetric  with  respect  to  the  z  axis,  let's 
consider  u  along  the  x-axis  which,  in  view  of  the  first  of  Eq .  (3),  permits 
the  expression: 

,  ,  2  2,1/2  .  x2  3d  -  v)t  S  2 

uj(x,0)  =  r  -  (r  -  x)  =  ~-  =  — - -  x  (4) 

Ed 

By  use  of  Eq.  (4),  the  film  stress  S  can  be  evaluated  from: 

1.  The  end  deflection  6c  of  a  cantilever  beam  of  length  t 

2.  The  relative  deflection  6  between  the  center  and  the  support  of  a 
simply  supported  beam  of  length  it 

3.  The  relative  deflection  6  .  between  the  center  and  the  support  of  a 
simply  supported  disk  of  diameter  i 

In  particular,  for  a  cantilever  strip  the  average  film  stress  is 


30  -  v)ti£ 


Further,  6^,  and  <$sd  for  the  same  film  stress  are  related  by 


6  =46  =  46 

c  s  sd 


wm 


/ 

V 


I 


fi 


t' 

\r‘ 

v 

i- 

1 


/ 


Next  consider  a  beam  clamped  at  both  ends  and  a  circular  disk  clamped 
along  the  periphery.  Since  the  substrate  is  subjected  to  isotropic  bending, 
the  application  of  additional  moment  at  the  boundaries  to  satisfy  the  zero- 
slope  condition  makes  the  deflection  of  the  entire  beam  or  disk  vanish. 


The  intrinsic 

stress 

distribution  across 

the  film  thickness  can 

be 

1 

determined  by  use 

of  the 

above  equations. 

Suppose  that  the 

intrinsic 

stress  o„  in  the  film  depends  on  the  distance  t 

from  the  substrate. 

Then 

the 

.  * 

average  stress  S  is 

related 

to  of  by 

tS(t)  =  /  o^(s)ds 


(7) 


If  at  a  given  point  does  not  vary  with  the  film  thickness,  we  can  differen¬ 
tiate  Eq.  (7)  to  obtain 


of(t)  =  S(t)  ♦  tS'(t) 


(8) 


The  expression  (Eq.  8)  enables  us  to  find  the  distribution  of 
of  approximately  if  the  average  stress  S(t)  is  measured  for  several  film 


thicknesses. 


Stresses  in  the  film-substrate  assembly  consistent  with  Stoney's  formula 
are  written  from  Eqs .  (1),  (2),  and  (3)  as  follows. 


Film  (superscript  f) 


(z)  =  o*'  (z)  =  of'(z)  S 
xx  yy 


"ij  ■ 0 


otherwise 


(I  <  z  <  |  +  t) 


Substrate  (superscript  S) 


(9) 


s  ,  .  s  ,  >  s,  . 

o  (z)  =  o  (z)  =  o  (z) 
xx'  yy 


6St  /  d i 

T*  lz  +  aJ 


°ij  =  ° 


otherwise 


(-  |  <  z  <  |) 
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Based  on  the  assumptions  and  Eqs .  (1)  through  (3),  one  can  also  obtain 
the  stress  field  at  the  f lim-substrate  interface  which  is  valid  at  the 
interior  of  the  assembly. 
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(10) 


However,  at  the  free  end  of  the  substrate,  there  is  a  boundary  layer  on  the 
order  of  the  film  thickness  in  extent,  within  which  non-zero  normal  and  shear 
stresses  exist  over  the  film-substrate  interface.  The  exact  solution  to  this 
problem  is  obtainable  by  use  of  classical  theory  of  elasticity  though  the 
analysis  is  expected  to  be  involved  and  time-consuming. 

Reference  5  provides  the  behavior  of  the  normal  and  shear  stress  fields 
within  the  region  of  the  boundary  layer  for  a  single  lap  joint  under  the 
application  of  uniaxial  tension  at  far  field.  Although  the  exact  stress 
distributions  at  the  end  of  the  film-substrate  interface  within  the  boundary 
layer  differ  from  that  in  a  single  lap  joint,  the  general  singularity  is 
similar.  In  both  cases,  the  normal  bond  stress  for  the  lap  Joint  is  highly 
localized  and  changes  sign  with  its  net  force  over  the  entire  bond  line 
vanishing  as  it  should.  A  similar  normal  stress  distribution  is  expected  for 
the  current  problem.  Shear  stress  on  the  interface  is  also  expected  to  be 
highly  localized. 

The  above  derivation  is  valid  when  the  thickness  of  the  substrate  is 
small  enough  such  that  it  will  undergo  a  change  of  radius  of  curvature.  When 
the  substrate  is  thick  enough  such  that  no  measurable  substrate  deformation 
can  be  made,  there  is  a  concern  as  to  the  magnitude  of  the  intrinsic  stress  in 
the  coating  in  comparison  to  that  in  a  thin  substrate. 

At  z  =  d/2,  the  plane  where  the  film  is  deposited,  the  strain  e(z)  is 
expressed  as 
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This  strain  needs  to  be  removed  when  the  substrate  is  very  thick.  Conse¬ 
quently,  it  is  equivalent  to  apply  a  stress  of  AS^.  of  a  magnitude  equal  to 

AS  =  |  J  (12) 

1  3  1  -  v  r 

When  this  compares  with  the  average  film  stress  S,  it  is  found  that 

4Sf  =  4  "3) 

which  implies  that 


ASf  =  o  (£  S) 


t  f  E 

since  -  <<  1,  unless  ^ - —  is  substantially  larger  than  - 
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The  small 


"o"  symbol  is  used  to  indicate  that  the  ratio  of  ASf  to  -S  is  very  small. 
Therefore,  it  is  concluded  that  the  dependence  of  the  intrinsic  stress  upon 
the  substrate  thickness  is  negligible  unless  the  film  modulus  is  much  greater 
than  the  substrate  modulus. 

B.  INTRINSIC  AND  THERMAL  STRESSES  IN  SUBSTRATE  DURING  COATING  DEPOSITION 


Both  intrinsic  and  thermal  stresses  develop  during  a  film  deposition. 
The  intrinsic  stress  is  due  to  the  atomic  interaction  between  the  substrate 
and  coating  materials.  The  thermal  stress  develops  due  to  the  temperature 
gradient  associated  with  the  deposition  process.  The  temperature  difference 
across  the  substrate  is  generally  very  small  (<  0.1°C)  because  of  the  thinness 
of  the  substrate.  This  difference  is  generally  disregarded  by  assuming  a 
constant  temperature  throughout  the  substrate/coating  assembly  in  most  of  the 
analytical  work.  This  assumption  may  lead  to  an  incorrect  deformation 
prediction  for  the  intrinsic  stress  calculation.  This  is  because  the  in-situ 
or  post  deposition  measurements  include  both  the  intrinsic  stress  and  the 
thermal  stress  deformations.  If  the  latter  part  of  the  deformation  is  not 


negligible  relative  to  the  former  one,  an  inaccurate  intrinsic  stress  would  be 
predicted.  A  thermal  deflection  analysis  for  a  thin  film  deposited  on  a  plate 
under  a  uniform  heat  flux  was  performed  since  a  direct  prediction  of  the 
intrinsic  stress  induced  deformation  is  difficult  to  achieve.  Blackburn  and 
Campbell published  an  in-situ  displacement  measurement  of  LiF  deposited  on 
a  glass  cantilever  strip.  There  appeared  to  be  a  sudden  jump  at  the  beginning 
and  a  sudden  drop  at  the  end  of  deposition.  This  deviation  from  a  smooth 
curve  has  not  been  resolved.  It  was  hoped  that  the  discrepancy  could  be 
explained  analytically  so  that  this  data  could  be  used  to  correlate  with  our 
analytical  predictions. 

Temperature  distributions  in  the  substrate  during  and  after  coating 
deposition  are  well  represented  by  solutions  to  the  one-dimensional  heat 
conduction  problem  of  a  slab  bounded  by  two  planes,  z  =  0  and  z  =  d 
(Fig.  1).  At  the  initial  instant,  the  temperature  T(z,t),  where  t  is  the  time 
variable,  is  equal  to  the  ambient  temperature  TQ.  During  the  film  deposition, 
we  assume  that  a  constant  rate  of  heat  flux  qQ  per  unit  area  is  supplied  at 
the  front  surface  and  heat  radiates  at  the  rear  surface  into  the  surrounding 
medium  of  temperature  Tq  according  to  Newton's  law,  with  coefficient  of 
surface  heat  transfer  H.  Once  deposition  is  ended,  the  heat  is  assumed  to 
radiate  from  the  front  and  rear  planes  with  a  linear  initial  temperature 
distribution . 

To  simplify  the  formulation  of  the  analysis,  we  introduce  the  following 
dimensionless  quantities: 

S  =  z/d 

t  =  <t/d2 

e  =  dH/K 

vU,t)  =  [T( z, t)  -  Tq]/[T1  -  T0] 

Bq  =  dq0/[K(T1  -  Tq)]  (IK) 
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where  K  is  the  substrate  conductivity,  <  is  the  substrate  diffusivity,  and  T1 
denotes  the  equilibrium  temperature  at  the  front  surface  when  the  deposition 
is  continued  indefinitely. 

In  this  dimensionless  formulation,  the  heat  conduction  equation 


3  v  3v 
_ .  2  =  3t 


(0  v  5  <  1, 
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is  accompanied,  for  the  deposition  problem,  by  the  boundary  conditions 
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and  similarly,  for  the  post-deposition  problem,  by 


v  (5,0)  =  a  +  b5 
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where  and  Vp  correspond  to  the  normalized  temperature  fields  for  the  depo¬ 
sition  and  post-deposition  problems.  For  the  detailed  derivation  of  solutions 
to  Eqs.  (15)  through  (17),  refer  to  Appendix  A. 

The  non-dimensional  temperature  fields  of  the  heated  surface  during  and 
after  deposition  are  depicted  in  Fig.  2  based  on  the  thermal  conductivity, 
specific  heat  of  the  substrate,  and  the  geometric  dimensions  of  the  cantilever 
strip.  The  temperature  distributions  versus  substrate  thickness  are  depicted 
in  Fig.  3  for  several  different  time  values.  They  correspond  to  a  temperature 
difference  between  the  heated  and  the  unheated  surfaces  of  0.0061°C  at  t  =  0.1 
sec  and  0.0096°C  at  t  =  100  sec. 


The  thermal  deflection  can  be  predicted  once  the  temperature  history  in  a 

substrate  is  determined.  To  calculate  thermal  deflection,  it  is  assumed  that 

a  circular  wafer  of  radius  a,  thickness  d,  and  coefficient  of  thermal 

expansion  a,  experiences  a  temperature  variation  T(z)  across  the  thickness 

over  a  concentric  circular  region  of  radius  b.  It  is  also  assumed  that  the 

t  7  \ 

boundary  is  simply  supported.  From  the  known  solution,'17  the  deflection  at 
the  center  is  obtained  as 

Ah2  d 

w(0)  =  [i  -  (1  +  v)  log(b'a)]  J  T(z)(z  -  d/2)  dz .  (18) 

dd  o 

The  tip  deflection,  6,  of  a  cantilever  plate  of  length  l  with  temperature 
distribution  T(z)  across  the  thickness,  is  known  to  be 

,  .2  d 

6  r  -  f  T(z)(z  -  d/2)  dz.  (19) 

d5  o 

The  integral  in  Eqs.  (18)  and  (19)  can  be  written  in  terms  of  the  dimension¬ 
less  temperature  (Eq.  14)  as 

d  -  1 

J  T(z)(z  -  d/2)  dz  =  (T.  -  Tn)  dr  /  vU,t)U  -  1  )/2)  d£ .  (20) 

0  u  0 

Thermal  deflection  based  on  Eq.  20  is  shown  in  Fig.  4;  also  plotted  are 
the  experimental  curve  for  the  total  deflection  by  Blackburn  and  Campbell  and 
the  difference,  which  corresponds  to  the  deflection  due  to  the  intrinsic 
stress.  The  thermal  deflection  curve  exhibits  a  sharp  initial  rise  and  drop 
at  the  end  of  the  deposition  period;  the  corresponding  sharp  changes  in  the 
total  deflection  curve  are  displayed  in  the  inset  with  a  magnified  scale. 

The  intrinsic  stress-induced  deflection  curve  is  smooth  and  free  from 
sudden  rises  and  drops.  It  also  places  the  intrinsic  stress  in  tension  over 
the  entire  time  domain.  The  calculated  thermal  deflection,  based  on  the 
temperature  gradient  of  a  constant  heat  flux,  predicts  a  rapid  rise  at  a  very 
short  time  after  the  deposition  begins.  The  calculation  explains  the  strange 
deflection  behavior  measured  and  reported  by  Blackburn  and  Campbell.  A  smooth 
intrinsic-stress-induced  deflection-history  curve  can  be  obtained  by  applying 


our  predicted  thermal  deflection  curve.  The  final  curve  not  only  permits  the 
deflection  to  approach  asymptotically  to  zero  at  small  time,  but  also 
eliminates  the  sudden  deflection  drop  at  the  end  of  deposition.  The 
analytical  prediction,  together  with  the  observation,  strongly  suggests  that 
the  small  temperature  difference  and  yet  large  thermal  gradient  between  the 
heated  and  unheated  surfaces  cannot  be  neglected  in  the  deformation  calcula¬ 
tion  for  the  intrinsic  stress  determination. 


III.  DEVELOPMENT  OF  DEFORMATION  MEASURING  TECHNIQUES 


The  accurate  prediction  of  the  intrinsic  stress  in  a  thin  film  relies 
heavily  on  the  ability  to  measure  the  deformation  of  the  system.  Therefore,  a 
viable  means  for  in-situ  and  post-vacuum  measurements  of  the  strains  must  be 
developed.  Furthermore,  the  technique  is  also  needed  in  the  determination  of 
the  material  properties  which  will  be  discussed  later  in  the  section. 

The  candidate  approaches  to  measuring  the  displacement  of  the  coating 
substrate  were  identified  as: 


1.  Laser  interferometry  (in-situ,  post-vacuum) 

2.  X-Ray  Diffraction  (post-vacuum) 

3.  Differential  capacitance  (in-situ,  post-vacuum) 

A.  LASER  INTERFEROMETRY 

The  Michelson  Interferometer  (see  Fig.  5)  consists  of  a  laser  source 
emitting  coherent  light  (indicated  by  parallel  lines  in  Fig.  5),  a  beam 
splitter,  a  reference  flat  reflector,  a  sample  substrate,  and  a  photo  plate. 
Light  from  two  separate  paths  both  impinge  on  the  photo  plate  and  create  an 
interference  pattern.  The  first  path  is  from  source  to  beam  splitter,  where 
it  is  reflected  to  the  reference  flat  and  then  back  through  the  beam  splitter 
to  the  photo  plate.  The  second  path  is  from  the  source  through  the  beam 
splitter  to  the  sample,  where  it  is  reflected  back  to  the  beam  splitter  and 
then  reflected  to  the  photo  plate.  The  nature  of  the  interference  pattern  is 


dependent  on  the  difference  in  path  length  between  the  two  optical  paths.  The 
photo  plate  is  always  illuminated  uniformly  as  long  as  the  sample  is 
absolutely  perpendicular  to  the  optical  axis  and  the  difference  between  the 
two  path  lengths  is  an  integral  multiple  of  the  wavelength  of  the  laser  light 
(Fig.  5).  If  the  path  length  at  the  sample  is  changed  by  one-half  wavelength 
by  moving  the  sample  along  the  optical  axis,  the  photo  plate  will  become 
uniformly  dark  because  the  two  light  waves  will  be  out  of  phase  and  will 
cancel  each  other.  As  the  sample  is  moved  along  the  optical  axis,  the  photo 
plate  will  be  alternately  uniformly  bright  or  uniformly  dark  depending  on 


whether  the  two  waves  reinforce  or  cancel  each  other. 
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Figure  6  illustrates  the  case  where  a  flat  sample  is  tilted.  Dark 
fringes  appear  at  the  photo  plate  where  the  peaks  (solid  lines)  interact  witn 
the  troughs  (dotted  lines)  as  the  sample  beam  interferes  with  the  reference 
beam.  The  separation  between  two  bright  or  dark  fringes  represents  one  half 
wavelength  of  displacement  at  the  sample.  (A  one-half  wavelength  displacement 
causes  the  optical  path  length  to  change  by  a  full  wavelength.) 

The  fringes  appear  as  concentric  circles  when  the  sample  curvature  is 
spherical  (concave  or  convex)  (Fig.  7).  To  test  whether  the  sample  is  concave 
or  convex  it  can  be  moved  slightly  towards  the  source.  If  the  circles  expand 
then  the  sample  is  convex.  If  they  contract  it  is  concave. 


When  a  thin  film  is  deposited  on  a  circular  substrate,  the  resulting 
stress  usually  causes  the  substrate  to  bend  spherically,  giving  a  pattern  of 
concentric  rings.  The  important  parameter  in  this  case  is  the  displacement, 
or  deflection,  at  the  center  relative  to  the  edge.  The  displacement  is 
measured  along  the  central  axis  perpendicular  to  the  substrate  and  is  the 
distance  the  substrate  surface  moved  as  a  result  of  stress  in  comparison  to 
its  pre-stressed  condition.  This  distance  can  be  measured  at  the  photo  plate 
by  counting  the  fringes  from  the  center  of  the  substrate  to  the  edge.  The 
displacement  is  equal  to  the  number  of  fringes  multiplied  by  one  half  the 
wavelength  of  the  laser  light.  In  actual  practice,  it  is  difficult  to  mount 
the  substrate  exactly  perpendicular  to  the  optical  axis.  In  this  case,  where 
the  substrate  is  at  a  slight  angle,  the  fringe  count  from  the  center  to  one 
edge  will  be  different  than  the  count  from  the  center  to  the  opposite  edge.  A 
simple  correction  factor  is  based  on  the  geometric  construction  shown  in 
Fig.  8.  The  corrected  displacement  is: 


6  =  F,  - 

c  1 


D1(F1  ~  F2> 

D1  +  D2 


(21) 


where 

F i  is  the  fringe  count  from  center  of  the  ring  system  to  one  edge  of 
the  coated  substrate. 

is  the  fringe  count  from  the  center  of  the  ring  system  to  the 
opposite  edge. 
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D1  is  the  distance  from  the  center  of  the  ring  system  to  the  same  edge 
used  for  F 1 . 

D2  is  the  distance  from  the  center  of  the  ring  system  to  the  edge  for 
F2 ' 

A  small  part  of  the  outer  rim  of  the  substrate  is  uncoated  due  to  the 
masking  effect  of  the  substrate  holder  during  film  deposition.  To  correct  for 

/  O  \ 

this  condition  the  following  equation  is  used'  ' 

r  =  r'(  1  -  AD/D)  ( AD  <<  D)  (22) 


where 

r'  is  the  radius  of  curvature. 

D  is  the  diameter  of  the  substrate. 

AD  is  the  difference  between  the  substrate  diameter  and  the  film 
diameter . 

It  is  this  corrected  value  for  r  that  is  used  in  Eq .  (3)  to  calculate  the 

stress . 

Additional  corrections  must  be  used  in  cases  where  the  film  is  thick  in 
comparison  to  the  substrate  (t  >  5)1  d).  These  and  other  considerations  are 
discussed  in  Ref.  4. 

Initially,  a  circular  wafer  was  selected  as  the  substrate  configura¬ 

tion.  A  substrate  mount  was  then  designed  to  accommodate  a  two-inch  diameter 
substrate  for  both  deposition  and  for  pre-  and  post-deposition  assessments. 
The  mount  was  designed  to  fit  in  the  sputter  deposition  chamber,  the  laser 
interferometer,  and  in  an  x-ray  diffraction  instrument  which  was  used  to 

obtain  additional  stress  assessment  data  for  correlation. 

A  general  layout  diagram  of  the  substrate  mount  is  shown  in  Fig.  9.  The 
substrate  is  shown  here  in  a  tightly  clamped  condition  with  a  differential 
fluid  pressure  causing  it  to  bow  outwards.  The  ability  to  exert  fluid 

pressure  on  the  substrate  was  used  to  assess  the  interferometer 's  ability  to 
measure  a  small  deflection  on  the  substrate.  Provision  was  also  made  to  study 
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Fig.  9.  Substrate  Mount 


the  feasibility  of  having  a  capacitance  bridge  to  measure  the  substrate 
deflection  by  using  a  rigid  conducting  plate  separated  from  the  substrate  by  a 
dielectric  spacer.  Electrical  contact  to  the  substrate  was  accomplished  by 
depositing  an  aluminum  film  on  the  upper  surface  of  the  dielectric  spacer. 
The  capacitance  measurement  technique  will  be  discussed  in  the  next  section. 

An  analysis  was  performed  to  predict  the  normal  differential  pressure 
versus  the  maximum  deflection  relation  for  a  thin  circular  silicon  wafer  with 
a  clamped  edge.  This  was  done  for  the  purpose  of  calibrating  the  laser 
interferometric  measurements.  Based  on  the  design,  a  5.08  cm  (2  in.)  diameter 
silicon  wafer  is  clamped  between  two  stainless  steel  plates  with  a  0.238  cm 
(3/32  in.)  wide  washer  along  its  edge.  The  central  part  of  the  wafer  has  a 
diameter  of  2.302  cm  (0.90625  in.)  and  is  free  to  deform  under  both  the 
external  load  and  residual  stress  from  the  coating. 

It  is  known  from  the  strength  of  materials  that  the  deflection  in  a  flat 
disk  under  the  action  of  external  load,  e.g.,  residual  coating  stress  or 
external  pressure,  depends  strongly  on  the  disk  thickness  and  the  wafer 
boundary  conditions.  When  the  wafer  is  circumferentially  clamped  and  is 
thick,  the  bending  stress  dominates  in  the  deflection.  The  membrane  stress 
becomes  more  important  however,  as  the  wafer  thickness  is  reduced.  Therefore, 
the  objective  of  the  analysis  was  to  investigate  the  deflection  sensitivity, 
and  the  contribution  from  membrane  and  bending  stresses  of  a  clamped  wafer  for 
various  thicknesses  and  pressure  levels. 

The  equation  for  determining  the  pressure  exerted  on  a  circular  wafer 
with  a  clamped  edge  as  a  function  of  maximum  wafer  deflection  is  expressed 


where 


P 


64 x  e^c-i  8  E  t 

3  lR  ‘  +  3  O  -  v)  R  J 


-  normal  differential  pressure  across  the  thickness  of  the 
wafer 


(23) 


R 


=  radius  of  the  wafer  to  the  clamped  edge 


J 


Young' s  modulus 


Poisson's  ratio 


plate  stiffness,  defined  as  Et^/12(1  -  v^) 


The  first  term  on  the  right  hand  side  of  the  equation  represents  the 
pressure  carried  by  bending  action,  while  the  second  term  represents  the 
pressure  carried  by  membrane  action.  It  can  be  seen  from  the  above  equation 
that  the  membrane  term  becomes  large  when  large  deflections, 
i.e.,  6^ /R  ratios,  prevail.  Since  it  is  our  goal  to  ensure  that  the 
deflection  of  coated  wafers  can  be  measured  with  good  sensitivity,  it  would  be 
advantageous  to  use  thin  wafers  as  long  as  they  do  not  undergo  undesirable 
distortion  due  to  handling. 

The  pressure  was  determined  as  a  function  of  the  5c/R  ratio  for  0.020  and 
0.038  cm  (0.008  and  0.015  in.)  thick  wafers  based  on  a  value  of  180  GPa  for 
E/(1  -  v)  and  a  Poisson's  ratio  of  0.4  for  silicon.  The  radius  was  based  on 
the  2.302  cm  (0. 90625  in.)  diameter  value. 

Figure  10  shows  that  the  bending  effect  is  dominant  with  a  6c/R  ratio  of 
0.01  or  less  for  a  0.038  cm  (0.015  in.)  thick  wafer.  For  a  6  /R  ratio  greater 
than  0.03,  the  membrane  part  controls.  They  are  both  significant 
for  «c/R  ratios  between  0.01  and  0.03.  In  the  case  of  a  0.02  cm  (0.008)  in. 
thick  wafer,  the  corresponding  limiting  Sc/R  values  reduce  to  0.005  and  0.02, 
respecti vely .  Also  a  thick  (0.038  cm)  wafer  requires  a  much  higher  pressure 
to  generate  the  same  deflection  as  compared  to  a  thin  (0.020  cm)  wafer.  Thus, 
higher  measuring  sensitivity  is  achieved  with  thin  wafers  for  a  given 
measuring  capability. 

The  original  plan  was  to  clamp  a  substrate  in  the  substrate  mount  and 
then  perform  deposition  and  stress  analysis  without  removing  the  substrate 
from  the  mount.  In  this  configuration,  the  thin  film  must  have  a  smaller 
radius  than  tue  substrate  because  the  bending  moment  takes  place  at  the  edge 
of  the  film  and  there  would  be  no  deflection  if  the  edge  of  the  film  was  at 
the  clamped  edge  of  the  substrate.  The  equation  for  this  configuration  is: 

4  X  ^c 
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Fig.  10.  Normalized  Pressure  Contribution 
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R  ^  radius  of  the  substrate  to  the  clamp 

rQ  -  radius  of  the  film 

d  =  substrate  thickness 

t  =  film  thickness 

The  optimum  film  radius,  rQ,  occurs  when  r0/R  =  1//e  where  e  is  the  natural 
number  2.728,  or  when  the  film  radius  is  about  .6  of  the  substrate  radius. 
This  method  is  less  sensitive  than  that  of  Eq.  (3)  by  a  factor  of  about  6.5, 
even  when  the  optimum  film  radius  is  used. 

The  deflection  vs.  pressure  of  a  silicon  substrate  with  a  5.08  cm  (2  in.) 
dia. ,  and  0.028  cm  (0.011  in.)  thickness  was  measured  by  clamping  the 
substrate  in  the  mount  and  varying  the  applied  air  pressure.  The  air  pressure 
in  the  mount  was  varied  by  means  of  a  bellows  type  pippette  pump  and  monitored 
with  a  glass  U-tube  manometer  containing  water.  Deflection  (from  the 
interferograms)  was  calculated  by  using  Eq.  (22).  The  data  are  as  follows: 


Pressure 

( inches  of  water) 

Deflection 

(pm) 

Difference 
Zero  Pressure 

from 
( pm) 

+  1  .0 

7.78 

-5.56 

0.0 

13.34 

-1 .  u 

17.95 

+4.61 

The  average  displacement,  then,  for  a  pressure  change  of  one  inch  of  water  was 
5.08  urn  with  an  uncertainty  of  i  .5  pm  or  ±  10?.  Note  that  at  zero  pressure, 
there  was  a  displacement  of  1 3 . 3^  um.  This  displacement  was  primarily  due  to 
pressure  exerted  by  the  clamping  ring  of  the  mount. 

The  expected  theoretical  center  displacement  is  found  to  be  4.5  pm  for  a 
5.08  (2  in.)  diameter  silicon  wafer,  0.028  cm  (11  mils)  thick,  (E  =  11.37  GPa 


and  \) 


r vw«" v* 


-  .A).  The  first  experimental  result,  at  2.54  cm  (  +  1.0  in.)  of  water, 
was  about  2 '%  above  the  theoretical  value,  while  the  result  at  -2.54  cm 
(-1.0  in.!  of  water  was  close  to  theoretical.  The  reason  for  a  higher  value 
at  positive  pressure  is  not  known  at  this  time.  An  uncertainty  of  t  .5  ym 
represents  a  stress  variation  of  ±  4.12  x  10  dyne/cm  for  a  0.028  cm  thick, 
5.08  cm  diameter  silicon  wafer  with  a  1  ym  thick  deposited  film. 

Two  thin  film  coated  substrates  were  evaluated  using  the  interferometer 
method.  One  consisted  of  Ge  on  a  Si  wafer  and  the  other  was  on  Si. 
Both  consisted  of  1  urn  thick  film  on  0.028  cm  (0.011  in.)  wafers.  Results 
were  as  follows: 


DEFLECTION  OF  SUBSTRATE 

Bef  ore 

After 

Fi  1m  ~ype 

Deposition 

(ym) 

Deposition 
(  ym) 

Dif ference 
(ym) 

Stress  (Eq  3 
(dyne/crrr ) 

Ge 

1.58 

53.95 

52.37 

4.40  x  109 

A1203 

4.46 

20.80 

16.34 

1 .42  x  109 

The  direction  of  the  stress  was  compressive  in  both  cases.  After  deposition, 
the  Ge  film  was  etched  away  in  a  sulfuric  acid  bath.  The  deflection  of  the 
substrate  was  then  2.06  urn.  The  value  agrees  with  the  measurement  prior  to 
coating  and  the  deflection  measured  after  deposition  was  due  primarily  to  thin 
film  stress. 

The  interferograms  made  before  and  after  coating  the  substrate  are  shown 
in  Fig.  11.  As  can  be  seen,  the  substrates  in  their  original  form  are  not 
flat  but  have  irregular  surface  features.  The  deviation  from  flatness, 
however,  is  small  in  comparison  to  the  deviation  after  the  coating  process. 
The  fringes  after  coating  are  so  close  together  towards  the  edge  of  the 
pattern  that  they  merge  and  cannot  be  counted  on  the  print.  The  original 
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Fig.  11.  Interferograms 


regat  i  v-.'  .  3  sharp  enough  :  ;•  resolve  t  n*  ;  r  .  ng'-s  right  tc  me  ej«n 

patter;:  and  car.  be  easily  cour ted  with  tue  a ic  of  a  microscope. 

An  unco ateci  substrate  was  measured  before  and  after  heat  tr  •  3 


determine  if  there  was  hysteresis  in  the  amount  jf  deflection.  The  deflection 
at  the  center-  of  the  substrate  before  heat  treatment  was  found,  by  laser 
interferometry ,  to  be  22.33  cm.  The  substrate  »as  then  placed  in  a  dO'-C 
oven  for  2.5  hours  and  then  cooled  to  room  temperature.  The  deflection  was 
then  measured  to  be  22.68  urn,  indicating  negligible  temperature  hysteresis. 
The  difference,  0.35  urn,  is  within  the  experimental  accuracy  of  the  inter¬ 
ferometer  measurement. 


An  uncoated,  but  clamped,  substrate  was  tested  in  the  interferometer  in  a 
stainless  steel  mount  with  an  attached  thermocouple  in  order  to  assess  the 
temperature  stability  of  the  system.  While  periodically  recording  the  fringe 
pattern,  the  mount  was  heated  by  radiation  from  a  heat  lamp  and  then  cooled 
back  to  room  temperature.  No  fluid  pressure  was  used.  The  displacement  as  a 
function  of  temperature  is  plotted  in  Fig.  12.  The  result  indicates  that  the 
displacement  changes  with  temperature  and  that  this  would  be  a  problem  for  in- 
situ  stress  evaluations  which  are  planned  as  a  later  phase  in  this  project. 
The  changes  in  displacement  are  due  to  a  GTE  mismatch  between  the  substrate 
and  its  mount.  For  this  reason,  a  new  substrate  mount  was  constructed  of 
Stackpol  e  230 1  GTE  graphite  which  matches  the  expansion  coefficient,  of  silicon 
very  e.oseiy. 

B  •  X - : ■  A'i  DIFFRACTION 

X-ray  diffract icr.  (XRD)  methods  have  commonly  been  used  in  the  measure¬ 
ment  of  stress  in  crystalline  materials. * Many  of  these  techniques 
determine  stress  by  measuring  a  lattice  parameter  of  a  material  in  several 
directions.  Other  XRD  methods  measure  stress  in  thin  films  by  determining  the 

radius  of  curvature  tnat  a  thin  film  imposes  on  a  single  crystal 
(  1  1  _  m  \ 

substrate.  Most  of  these  techniques  employ  the  use  of  a  Lang  camera  or 
other  special  instrumentation;  however,  one  met hod ^  ^  ^  makes  use  of  a  simple 
version  of  a  double-crystal  diffractometer.  Its  use  for  making  radius  of 
curvature  measurements  is  simple,  provided  a  double-crystal  diffractometer  is 
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curvature  of  several  samples  were  analyzed  by  the  XR1- 


Tie t hod  and  were  also  determined  using  the  laser  interferometer.  This  served 
as  a  cheek  or.  the  two  methods  and  provided  additional  data  for  evaluation  in 
the  selection  of  a  stress  measurement  method. 

Radius  of  curvature  was  determined  using  a  modified  version  cf  the 
technique  cited  in  Ref.  11.  Measurements  were  performed  on  a  commercial ly 
available  (Blake  Industries)  double-crystal  diffractometer  using  copper  radia¬ 
tion  and  a  scintillation  detector.  The  diffractometer  was  configured  in  the 
non-parallel  unequispac ing  asymmetric  (  +  ,  -)  setting  (Fig.  13). A  (100) 
single  crystal  silicon  wafer  was  used  as  a  monochromator  and  was  aligned  so  as 
to  diffract  CuKa  radiation  from  the  (311)  lattice  planes  in  the  wafer.  The 
purpose  for  selecting  the  (311)  reflection  was  to  expand  the  x-ray  beam  so 
that  it  would  cover  as  large  an  area  as  possible  on  the  sample.  This  expan¬ 
sion  is  possible  because  one  set  of  (311)  planes  meets  the  (100)  surface  at  an 
angle  which  is  2.8°  less  than  the  Bragg  angle  (0).  In  this  configuration, 
the  incident  x-rays  strike  the  surface  of  the  monochromator  at  a  very  low 
angle  (2.8a),  in  order  for  diffraction  from  the  (311)  planes  to  occur,  and 
thereby  cover  a  large  area. 

Samples  consisted  of  0.028  cm  thick  5.08  cm  diameter  (100)  silicon  wafers 
which  were  aligned  to  diffract  from  the  (100)  planes  [(400)  reflection].  The 
sample  holder  of  the  diffractometer  was  modified  to  accept  the  specially 
constructed  substrate  holder.  During  initial  alignment  of  the  sample,  an 
approximately  3  mm  x  3  mm  aperture  (not  shown  in  Fig.  13)  was  placed  between 
the  monochromator  and  the  sample. 

Alignment  of  the  sample  consisted  of  monitoring  the  intensity  of  the 
diffracted  beam  while  the  sample  was  tilted  about  a  horizontal  axis  and 
rotated  ("rocked")  about  a  vertical  axis.  The  sample  tilt  was  adjusted  until 
it  produced  the  greatest  intensity  when  "rocked".  At  this  point,  the  normals 
of  the  diffracting  (311)  lattice  plane  in  the  monochromator  and  the  (400) 
lattice  planes  in  the  sample  are  both,  as  nearly  as  possible,  in  the  same 
plane.  This  will  reduce  the  rocking  curve  half  width  of  the  sample  to  the 


point  where  it  is  only  controlled  by  the  properties  of  the  material  ana  any 
external  strain.  Because  of  the  dispersive  geometry  used  in  the  diffrac¬ 


tometer,  both  the  Ka^  and  Ka^  reflections  will  be  observed  in  the  rocking 
curves. For  the  combination  of  Si  (311)  and  Si  (400),  the  Ka 1  -Ka^  peak 
separation  (dispersion  factor)  is  79-8  arc  seconds  for  Cu  radiation. 

Following  sample  alignment,  rocking  curves  were  recorded  with  two  1  .5  mm 
diameter  apertures  placed  between  the  monochromator  and  sample.  The  apertures 
and  samples  were  placed  so  that  the  two  x-ray  beams  struck  the  sample  at 
points  on  a  diameter  of  the  wafer  and  equidistant  from  the  center  (±  2  mm). 
An  unstressed  wafer  will  diffract  from  both  beams  ( A , B )  at  the  same  angular 
position  (Fig.  14).  In  a  stressed  wafer  both  beams  will  not  diffract  at  the 
same  angular  position  (Fig.  15).  The  separation  of  the  diffraction  peaks  from 
each  beam  is  then  a  measure  of  the  angular  difference  between  the  tangents  to 
the  points  the  two  beams  strike  the  sample.  A  series  of  rocking  curves  was 
run  on  each  sample  with  aperture  spacings  (L)  of  1.0, 1.4,  1.5  and  1.7  mm. 
Rocking  curves  were  run  with  different  aperture  spacings  as  a  check  to  see  if 
the  wafers  were  deformed  spherically  and  produced  the  same  radius  of  curvature 
independent  of  aperture  spacing.  In  one  instance,  rocking  curves  of  a  sub¬ 
strate  were  taken  before  and  after  film  deposition  in  order  to  calculate  the 
stress  produced  by  the  film.  The  sensitivity  of  the  technique  is  limited  by 
the  resolution  of  the  rocking  curve  and  the  distance  between  the  apertures. 

An  analysis  was  performed  to  determine  the  most  favorable  spacing  of  two 
parallel  x-rays  when  the  double  crystal  XRD  technique  is  used.  The  equation 
for  determining  the  slope  of  a  thin  circumferentially  clamped  wafer  within  the 
small  deflection  theory,  under  normal  pressure  is 


s  =  -3pn( 1  -  v2)(R2 


n2)/4Et3 


while  the  equation  for  the  slope  difference  of  two  points  at  a  distance 
apart  along  the  diametric  line  can  be  expressed  as 

ds  =  -3p( 1  -  v2)  L(3R2  +  3RL  +  L2  -  R2)/4Et3  (2 
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Fig.  14.  X-ray  Rocking  Curve  of  Unstressed  Silicon  Wafer 
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Fig.  15.  X-ray  Rocking  Curve  of  Sample  B,  L 


where 


s  =  slope 

p  =  uniformly  distributed  load  intensity 

R  =  radius  of  the  wafer 

n  =  radial  distance  from  the  center  of  the  wafer 

v  =  Poisson's  ratio 

E  =  Young's  modulus 

L  =  distance  between  the  two  fixed  x-rays 

For  this  analysis  the  ratio  of  s/p  vs.  r  was  calculated  to  eliminate  the 
need  for  a  special  load  intensity.  The  values  used  for  Young's  modulus  and 
Poisson's  ratio  were  that  of  silicon;  E  =  108  GPa  (15.66  msi)  and  v  =  0.4. 

The  slope  of  the  wafer  as  a  function  of  r  is  calculated  from  Eq.  (25).  A 
diametric  line  was  used  in  the  calculation  because  it  has  the  maximum  rate  of 
change  of  slope  of  any  chord  on  a  wafer.  The  slope  difference  was  obtained  by 
calculating  the  slopes  at  the  locations  of  rays  1  and  2  which  are  at  a  spacing 
L  apart  (Fig.  16)  and  then  taking  the  difference  of  the  two  slopes  Eq.  (26). 
The  procedure  was  performed  by  first  placing  two  imaginary  rays  symmetrical 
about  the  center  of  the  disk  and  calculating  the  slope  difference.  These 
imaginary  rays  were  then  traversed  relative  to  the  wafer  by  maintaining  the 
spacing  distance  L.  The  slope  difference  was  calculated  at  an  increment  of 
one-fiftieth  of  the  radius  until  ray  2  reached  the  disk  edge.  The  clamped 
wafer  was  assumed  to  have  an  inner  radius  of  2.302  cm  (0.90625  in.)  and  a 
thickness  of  0.028  cm  (0.011  in.)  for  the  analysis.  The  spacings  for  L  were 
1,  1.7  and  2  cm  (0.3937,  0.6693,  and  0.7874  in.). 

The  graphs  for  the  three  spacings  which  were  used  for  the  0.011  in.  thick 
wafer  are  shown  in  Figs.  17  through  19.  The  lower  curve  in  each  figure 
represents  the  slope  along  any  diametric  line  as  a  function  of  radial  distance 
from  center  of  the  wafer.  The  horizontal  coordinate  of  the  slope  difference 
curve  measures  the  location  of  ray  1  while  ray  2  is  at  a  positive  distance  L 
away.  In  all  three  cases  the  maximum  slope  difference  occurs  when  rays  1  and 
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Fig.  16.  Schematic  of  Wafer  Configuration 
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2  are  symmetrical  about  the  center  of  the  disk.  The  maximum  slope  differences 
are  0.0097,  0.0142,  and  0.017  for  the  three  cases  (1.0,,  1.7  and  2.0  cm) 

cons idered . 

The  radius  of  curvature  (r)  for  the  stress  calculation  in  Eq .  (3)  is 

obtained  from  the  relationship 


where 

A 9  =  the  amount  of  angular  rotation  between  the  reflections  from 

the  two  apertures 

L  =  the  distance  separating  the  two  x-ray  beams  on  the  sample 

The  direction  of  the  stress  can  be  determined  by  noting  whether  the  substrate 
is  concave  or  convex.  This  is  accomplished  by  alternately  blocking  one  of  the 
apertures  and  rocking  the  sample  through  the  x-ray  reflection.  If  the  sub¬ 
strate  is  convex  the  beam  A  (reference  Fig.  13)  will  be  recorded  first  as  e  is 
decreased.  The  opposite  is  true  if  the  substrate  is  concave. 

In  calculating  the  stress  in  a  thin  film,  the  laser  interferometer  method 
measures  the  displacement  (6c)  between  the  center  of  the  wafer  and  the  edge. 
To  allow  comparison  with  the  interferometer  method,  the  quantity  is 

calculated  from  the  radius  of  curvature,  as  determined  by  XRD,  with  the 

following  expression: 


«e-fc 


The  underlying  assumption  in  the  XRD  method  is  that  the  wafer  is 
originally  flat,  or  deformed  spherically,  and  that  it  deforms  spherically  as  a 
result  of  film  deposition.  This  is  because  the  XRD  method  only  measures  the 
angular  difference  between  two  points  and  then  calculates  R  and  d  based  on  a 
spherical  geometry.  It  will  be  shown  that  this  underlying  assumption  is  not 
always  true. 


•’’i 


Three  as-received  wafers  (A-C)  were  analyzed  by  both  the  XRD  and  the 
laser  interferometer  methods.  An  additional  wafer  (D)  which  had  a  strain 
gauge  mounted  to  its  back  surface  was  also  analyzed  by  both  techniques,  before 
and  after  deposition  of  a  sputter  deposited  Si  film.  Examples  of  the  XRD 
rocking  curves  are  given  in  Figs.  15,  20  and  21.  A  summary  of  the  XRD  and 
interferometer  data  is  given  in  Table  3-  The  data  indicates  that  the  as- 
received  silicon  wafers  are  not  flat,  or  that  the  sample  mount  may  be 
distorting  the  wafers,  since  deformations  (in  some  cases  highly  non-spher i cal ) 
are  present.  This  is  most  easily  seen  in  the  interferograms ,  where  concentric 
equally  spaced  circular  lines  would  be  evidence  for  spherical  deformation. 

In  some  instances  (samples  A,  B)  the  agreement  between  the  XRD  and  laser 
interferometer  data  is  reasonably  good  (A  <  10(£)  but  in  other  cases  (sample  D) 
the  differences  are  as  great  as  a  factor  of  two.  Several  factors  may  be 
responsible  for  the  discrepancies  between  the  XRD  and  interferometer  data. 
These  include: 

1.  The  XRD  method  assumes  that  the  wafer  is  spherically  deformed,  but 
interferograms  suggest  that  this  may  not  be  the  case. 

2.  The  locations  of  where  the  two  x-ray  beams  strike  the  sample  is 
assumed  to  be  exactly  known,  but  actually  is  only  known  to  +  2  mm 
(in  both  x  and  y).  This  would  affect  the  i nterpretat ion  of  the 
interferograms . 

3.  The  interferograms  are  also  sensitive  to  non-stress  related  height 
differences  (wafer  or  thin  film  thickness  variations)  that  are 
interpreted  as  being  a  result  of  stress. 

4 .  In  some  instances  the  Ka1  reflection  from  one  x-ray  beam  may  be 

partially  (Fig.  20)  or  totally  superimposed  (Fig.  21)  on 

the  Ka^  reflection  from  the  other.  This  may  produce  a  slight  peak 
shift  in  the  observed  position  of  the  Ka ^  reflection. 

The  fact  that  the  wafers  are  not  spherically  deformed  is  also  indicated  by  the 
XRD  data  since  the  calculated  radius  of  curvature  often  varies  considerably  on 
a  given  wafer  depending  on  the  separation  distance  (L)  of  the  two  x-ray  beams. 

The  sensitivity  of  the  XRD  method  is  limited  by  the  separation  distance 
(L)  of  the  two  x-ray  beams  on  the  sample  and  the  resolution  of  the  two  dif¬ 
fracted  x-ray  beams.  The  aperture  of  the  detector  limits  L  to  be  no  greater 


than  3-01  cm  and  the  resolution  of  the  rocking  curves  is  on  the  order 
of  5  arc  seconds.  The  largest  detectable  "adius  of  curvature  is  then  approxi¬ 
mately  1240  m.  For  a  film  thickness  (t)  of  0.10  ym  and  a  radius  of  curvature 
(r)  of  1240  m,  the  smallest  measurable  displacement  (6  )  is  280  pm  and  since 
for  {100}  Si : 

E 

— — - — r  =  1.805  x  101  dynes/cm^  (Ref.  13) 

(1  -  v3) 

Then,  from  Eq.  (1)  and  the  above  stated  conditions,  the  minimum  detectable 

Q  p 

thin  film  stress  S  =  1.89  x  10°  dynes /cm.  This  compares  with  a  value  of  6  x 
7  P 

10'  dynes/cm  reported  in  Ref.  14  using  Fe  radiation  and  the  (440) 
reflection.  It  is  to  be  noted  that  the  non-spher ical  deformation  places  a 
limitation  on  the  XRD  method  of  stress  measurement,  since  it  calculates  the 
radius  of  curvature  (r)  and  the  displacement  (6  )  based  on  a  spherical 
geometry.  A  further  limitation  with  the  double-crystal  XRD  method  of  stress 
measurement  is  its  dependence  on  single  crystal  substrates.  This  would 
eliminate  it  from  consideration  with  polycrystalline  or  amorphous  substrates 
which  are  planned  for  future  use  in  the  program. 

C.  DIFFERENTIAL  CAPACITANCE 

It  has  been  demonstrated  experimentally  that  laser  interferometry  can  be 
successfully  used  to  characterize  the  flatness  of  pre-deposi ted  wafers  and  the 
deformation  of  post-deposited  wafers.  The  sensitivity  of  the  displacement 
measurement  depends  upon  the  wavelength  of  the  laser.  However,  if  the  in-situ 
deformation  is  desired  during  the  coating  deposition,  problems  are  expected 
due  to  the  difficulty  of  setting  the  light  path  and  the  vibrational  motion  of 
the  sputtering  system.  It  was  determined  after  a  careful  study  that  laser 
interferometry  may  not  be  feasible  for  the  in-situ  deformation  measurement. 
The  double-crystal  XRD  technique  is  also  applicable  for  post  deposition 
deformation  measurement.  Like  the  laser  interferometer,  the  vibrational 
motion,  together  with  the  bulky  apparatus,  make  the  measurement  impractical. 
In  addition,  the  in-house  XRD  equipment  cannot  be  used  on  polycrystal  line  or 
amorphous  structures. 
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r'-  rent  :  a  I  capacitance  method  has  been  used  as  a  viable  way  to 

.deflect  .on  of  structures.  It  works  on  the  principle  that  the  elec- 

•  ■  o  •. .  eapae:  ranee  between  two  conductor  plates  separated  by  an  insulator  is 

linearly  inversely  proportional  to  the  distance  between  the  two  plates.  From 

past  o  a  •  ■  ce .  a  change  of  dimension  Ad  between  two  capacitance  plates, 

_c 

v v  arc  a:  distance  d  apart,  can  be  measured  to  10  d  by  the  use  of  the 


(29) 


wr.cr,  Z  is  it:  the  range  of  10  pF.  Therefore,  a  deflection  sensitivity  of 
0.  ’3  am  (5.0  uin.)  can  be  achieved  if  the  initial  gage  distance  between  the 
two  plates  is  no  more  than  0.013  cm  (0.005  in.).  The  normal  deflection  of  an 
edge-clamped  circular  wafer  under  the  application  of  normal  pressure  can  be 
measured  by  the  use  of  the  differential  capacitance  method  by  measuring  the 
capacitance  between  the  substrate  and  a  rigid  metal  plate.  The  differential 
change  in  capacitance  as  a  function  of  pressure  is  given  by  Neubert  as: 

AC  =  .0174  v  }  -  p  (30) 

Ed  t3 

w.te.’w  p,  K,  t  are  tne  applied  pressure,  radius  and  thickness  of  the  plates, 
respectively,  and  a  is  the  distance  between  the  plates.  The  initial 
capacitance  is  that  of  a  pair  of  circular  plates  separated  by  a  dielectric: 


where  A  : s  the  area  of  the  circular  plate.  The  fractional  change  of  capacity 

as  a  ft-nrO  , nr,  of  pressure  is  then: 


AC/C  -  0.0625 


P  4 
(i  -  / )  «r 


Combining  Eq.  (30)  with  Eq.  (23),  neglecting  large  deflection  effects,  we 


6c  -  3d  AC/C  (33) 

The  capacitance  between  the  Si  substrate  and  the  fixed  metal  plate  was 
measured  by  means  of  an  impedance  analyzer  (Hewlett  Packard  Type  4191),  The 
air  pressure  in  the  substrate  mount  was  controlled  and  monitored  as  it  was  for 
the  deflection  vs.  pressure  experiment.  A  test  frequency  of  4  MHz  was  used  in 
the  impedance  analyzer.  Capacity  resolution  was  about  ±  .2  pF.  The  total 
measured  capacity  was  made  up  of  the  following  components: 

1.  Substrate  initial  capacity  (inside  clamp)  117.7  pF 

2.  Capacity  of  clamp  ring  77.0 

3.  Capacity  of  connecting  cable  91.9 

Total  capacity  286.6  pF 

A  comparison  of  the  MSL  measured  capacitance  values  to  the  calculated  values 
based  on  Eq.  (3)  is  shown  below. 
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By  Eq.  (33?,  the  deflect ioi.  due  to  z  .2  pF  (the  approximate  uncertainty  of  the 
above  data)  :s  calculated  to  be  ±  .65  cm.  The  equivalent  stress  (assuming  a 
or.e  micron  deposited  film)  is  t  3.86  x  10°  dyne/enr  (5600  psi).  This  repre¬ 
sents  the  uncertainty  cf  the  capacitive  method  with  a  clamped  substrate.  The 
z  .  cB  „t  uncer-aintv  of  measured  deflection  is  not  much  worse  than  the  t 
.5  ..m  jr.eertair.ty  of  tne  interferometry  method  (Appendix  B,  page  B-1 ) ,  but  an 
arc. an. pod  substrate  can  DC  used  when  the  interferometry  method  is  used  for 
stress  measurements,  giving  an  order-  of  magnitude  improvement  of  stress 
s  t*  r  i  s  i  t  .  v  -  r.  y  . 

ln-situ  test  runs  were  made  to  measure  the  deflection  of  Si  wafer 
substrates  during  coating  deposition.  It  was  observed  that  RF  signal  inter¬ 
ference  was  rot  a  problem  during  the  entire  course  of  the  data  acquisition. 
Also,  the  vibrational  motion  of  the  sputtering  system  did  not  affect  the  mea¬ 
surement  since  the  capacitance  device  was  an  integral  part  of  the 
specimen  ''holder  assembly.  Consequently,  it  has  been  concluded  that  the  dif¬ 
ferential  capacitance  method  is  the  best  among  all  three  techniques  in  the 
measurement  of  substrate  deformation  for  pre,  in-situ,  and  post-deposition 
configurations.  Test  runs  were  made  for  bare  wafers  and  post-deposition 
wafers  with  4  am  thick  Si/Ge  discrete  and  rugate  coatings.  Table  4 
illustrates  the  comparison  between  the  laser  interferometry  and  the 
differential  capacitance  methods.  It  is  seen  that  the  results  of  using  two 
methods  differ  by  30  -  40£,  which  might  be  caused  by  possible  sample  movement 
as  a  result  of  sample  holder  handling. 

A  concern  was  raised  regarding  the  actual  boundary  (edge)  behavior  of  a 
circular  wafer  in  the  sample  holder.  It  has  been  predicted  analytically  that 
the  maximum  center  deflection  of  a  simply  supported  circular  plate  under 
normal  pressure  is  about  four  times  as  large  as  that  when  it  is  clamped  along 
its  edge.  Since  the  former  case  has  better  deflection  sensitivity,  the  edge 
was  only  "finger  tight"  in  the  actual  assembly  to  provide  the  best  possible 
simply  supported  boundary  condition  while  still  minimizing  possible  non- 
stressed  related  movement.  However,  the  "finger  tight"  boundary  condition 
will  not  eliminate  possible  specimen  movement. 
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It  was  found  during  the  calibration  runs  of  the  capacitance 
instrumentation  that  reasonable  data  were  obtained  when  the  wafer  was  placed 
horizontally  with  no  clamping  force  from  the  holder  and  with  only 
gravitational  force  exerting  downward.  However,  when  the  wafer  specimen  is 
placed  vertically  in  the  holder,  there  needs  to  be  some  clamping  force  (e.g., 
finger  tight)  from  the  holder.  The  magnitude  of  this  force  cannot  be 
measured.  Also,  there  might  be  a  small  gap  existing  between  the  wafer  and  the 
holder  in  some  places.  Consequently,  no  reliable  deflection  data  could  be 
obta ined . 


After  reviewing  the  other  supporting  methods  and  considering  that  high 
deformation  sensitivity  is  vitally  important,  it  was  decided  that  the  canti¬ 
lever  method,  rather  than  the  circular  wafer  method,  should  be  used  to  measure 
substrate  deflection.  The  decision  was  based  upon  the  following  observations: 


1.  The  tip-deflection  of  a  cantilever,  as  shown  by  Eq.  (6),  due  to 
film-stress  is  4  times  the  central  deflection  of  a  simply  supported 
circular  wafer  if  every  condition  is  the  same  and  the  length  of  the 
cantilever  is  equal  to  the  diameter  of  the  wafer. 


• 


con! :gur 


i'ne  specimens  used  for-  in-situ  measurement  need  to  be  placed  ver¬ 
tically.  The  cantilever  can  be  clamped  without  any  difficulty, 
while  the  simply  supported  condition  for  the  wafer  is  very  difficult 
to  achieve  in  the  vertical  position. 

“■ijc  of  the  published  data  on  film-stress  have  been  obtained  by  the 

cart i ! ever  method. 

Analytical  expressions  of  the  physical  quantities  are  much  simpler 
and  more  compact  for  the  cantilever  than  for  the  circular  wafer. 

Some  improvements  on  the  cantilever  method  have  been  proposed,  but 
improvements  to  the  circular  wafer  method  seem  to  be  difficult. 

ew  sample  holder  and  specimens  have  been  fabricated  for  the  cantilever 
at  ion.  More  calibration  runs  for  bare  and  thin-film  deposited  speci¬ 


mens  will  be  performed. 


_•  V  v  /  /  /  /  •• 


: V.  MATERIAL  PROPERTY  DETERMINATION 

>  the  assessment  of  the  stress  response  in  the  coating  for  a  thin- film 
deposited  suostrate,  the  material  properties  for  thin-film  materials,  e.g. , 
Young's  rod-Ius,  shear  modulus,  coefficient  of  thermal  expansion,  need  to  be 
determined.  Unfortunately,  the  properties  of  thin-film  materials  are  not 

C  1 7  \ 

necessarily  the  same  as  the  bulk  materials.  Also,  the  different  coating 
techniques,  e.g.,  evaporated,  electroplated,  and  chemically  deposited  coatings 
yield  different  properties.  Therefore,  an  accurate  prediction  of  the  film 
stresses  under  operational  environments  is  contingent  on  reliable  property 
input  for  the  coating  materials. 

Because  of  the  very  thin  character  of  thin-film  coatings  (hnoo  A  or 
less),  the  measurement  of  their  Young's  moduli  cannot  be  achieved  by  _ne  use 
of  conventional  methods.  Some  measurements  have  been  done  by  other  investiga¬ 
tors  using  the  bulging  of  a  coating  from  an  applied  normal  pressure  to  the 
coating  when  part  of  the  substrate  was  etched  away.  This  is  valid  only  if  the 
moduli  of  coatings  are  not  a  strong  function  of  the  intrinsic  stresses.  When 
the  properties  change  with  respect  to  the  existing  intrinsic  stress,  other 
techniques  need  to  be  explored. 

Cantilever  reeds  or  strips  are  being  used  for  the  substrate  configuration 
for  the  materia^  property  measurement.  The  coatings  will  be  deposited  with 
controlled  parameters  such  as  temperature,  deposition  rate  and  total  thick¬ 
ness  . 

A.  DETERMINATION  OF  YOUNG'S  MODULUS 


The  fundamental  bending  frequencies  of  the  specimen  before  and  after  the 
deposition  will  be  measured.  Due  to  the  film  deposition,  the  added  coating 
affects  the  fundamental  frequency  of  the  cantilever  substrate.  Since  the 
coating  is  an  integral  part  of  the  cantilever  substrate,  Young's  modulus  of 
the  coating  materials  car.  be  expressed  in  terms  of  frequency  change  and  other 
specimen  and  material  parameters. 
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I  t  is  realised  that  instr insic  stresses  develop  as  a  result  of  coat:  rig 
deposition.  The  magnitudes  of  the  intrinsic  stresses  may  be  large  enough  to 
cause  yielding  in  the  coating.  Therefore,  it  is  possible  that  the  tangent 
modulus  at  the  stress  level  differs  from  the  initial  tangent  modulus  (Young's 
modulus)  of  the  material.  One  may  thus  expect  that  the  measured  modulus  using 
a  vibration  cantilever  will  not  give  the  Young's  modulus  value.  In  fact,  when 
the  specimen  is  vitiated,  the  stress  in  the  coating  will  increase  first  to  its 
maximum  value  and  then  decrease  during  the  first  cycle.  Although  it  may 
follow  the  nonlinear  part  of  the  stress-stra in  curve  during  initial  loading, 
the  linear  unloading  stress-strain  curve  will  be  followed  in  the  subsequent 
cycles.  If  no  fatigue  damage  is  assumed  to  occur  from  the  vibration,  the 
Young's  modulus  of  the  coating  is  expected  to  be  obtained  from  the  vibration 
technique . 


The  basic  equation  for  the  fundamental  frequency  of  a  uniform  cross  sec¬ 
tioned  cantilever  strip  is  furnished  in  Ref.  1 8 , 


2  ,  _ 


where  E  p  ,  and  d  are  Young's  modulus,  density  and  thickness  of  the  sub¬ 


strate,  respectively,  l  is  the  unsupported  length  of  the  cantilever  and  aQ  is 


a  constant  (aQ  -  1.8751).  As  the  film  is  deposited  onto  the  substrate,  the 


fundamental  frequency  f  shifts  to  f  according  to 


-  f _ ^ _ ) 

2  ^  p  d  +  p  -d  „ ' 
l  s  s  f  f 


where 


A  =  |d3E  +  d3E  +  3drd  (d  +  dj2  E  Ep/(d  E  +  drE_}/12 

SS  II  ISS  I  SI  SS  II 


o,  -  density  of  the  film 
df  =  thickness  of  the  film 

=  Young's  modulus  of  the  film 


The  vibrating  reed  method  has  been  successfully  employed  to  determine  Young's 

/  iO\ 

moduli  of  thin  films  coated  on  substrates'  ’  and  to  study  the  adherent 

M  q  \ 

characters  of  the  coated  film.'  This  method  does  not  require  the  separa¬ 

tion  of  the  film  from  the  substrate  and  is  especially  useful  for  determining 
the  moduli  of  materials  that  are  stable  only  in  the  form  of  the  films  on  a 
substrate.  Rugate  films  containing  amorphous  Germanium  are  known  to  be  stable 
only  in  the  coated  form.^^) 

To  determine  Young's  modulus  of  the  deposited  coating,  a  good  under¬ 
standing  of  the  fundamental  frequency  change  versus  the  ratios  of  density  and 
thickness  of  the  two  materials  (substrate  and  coating)  is  important.  There¬ 
fore,  an  assessment  of  the  sensitivity  of  E^/Eg  vs.  p^/ps,  df/ds*  and  Af/fQ  is 
needed  for  the  successful  prediction  of  the  coating  modulus. 

Given  a  coating  thickness  such  that  df/ds  <<  1,  Eq.  (35)  reduces  to 


where 


d p  dp  2 

r  =  Mi  .  a  /  .  b  (/)  ) 


a  ■  1  tr>  -  5  £) 

s  s 


4  ^0  ;uE  ' 
s  s 


8 J T  8  - 
s  s 


Substitution  of  the  various  p^/ps  and  df/ds  values  into  Eq.  (36)  gives  an 
equation  of  Af/f  in  terms  of  Ep/Eg.  For  any  range  of  Ef/Eg  values,  the 


correspond ing  values  of  Af/f  can  be  found. 


A  substrate  of  Si  and  a  film  of  Ge  have  been  used.  Equation  (3o)  is  a 
general  formula  and  car.  be  used  for  any  combination  of  materials,  but  Si  arid 
Ge  are  currently  being  used  in  this  program.  Testing  has  been  ccnducted  on  Si 

wafers  0.02794  cm  (0.011  in.)  thick  with  a  density  of  2.33  g/cc.  Young's 

modulus  of  the  Ge  coating  is  the  element  of  interest  in  this  analysis;  its 
thickness  ranges  from  1,000  to  3,000  A  and  its  density  from  4.256  to  6.384 
g/cc  (see  Table  5  for  details).  The  coating  thickness  is  in  the  range  of 
existing  coating  thicknesses  while  the  20%  variation  is  used  to  cover  possible 
density  changes.  Three  thickness  ratios  of  3.58  x  10~\  7.16  x  10-i|,  and 

1.074  x  10~3,  and  five  density  ratios  of  1.83,  2.05,  2.28,  2.51,  and  2.74  are 
obtained  from  the  data  gathered.  When  these  values  are  substituted  into 
Eq .  (34),  one  obtains  fifteen  different  cases  for  the  Af/f  vs.  Er/E  rela- 
t ions'- ip.  The  graphs  of  the  fifteen  curves  are  plotted  as  three  separate 
graphs  with  the  density  ratio  varying  for  each  of  the  three  thickness 
ratios.  Each  of  the  three  graphs  has  been  plotted  with  Af/f  vs.  Ef/E„,  as 
depicted  in  Figs.  22,  23,  and  24.  The  uppermost  curve  represents  a  density 
ratio  of  1.83,  and  with  each  succeeding  curve,  the  values  are  2.05,  2.28, 

2.51,  and  2.74.  Figures  22,  23,  and  24  represent  thickness  ratios  of  2.38  x 
1 0_i4 ,  7.16  x  10”\  and  1.074  x  10"^,  respectively. 

Table  5.  Material  Properties 
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interpreted  from  the  equations  and  tne  thr 


grains: 


For  i,.d,  'v  1 ,  the  modulus  ratio  is  linearly  proportional  to  the 
density  ratio  and  the  frequency  ratio  ( A f  f  )  is  linearly  propor¬ 
tional  to  the  thickness  ratio.  For  large  d^-d  values,  the  higher' 
order  terms  of  dp/d  will  contribute  more. 
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B.  determination;  of  shear  modulus 

Two  independent  elastic  moduli  are  needed  to  characterize  the  elastic 
properties  of  an  isotropic  material.  In  this  respect,  it  is  essential  to 
characterize  the  mechanical  properties  of  a  coated  film  by  measuring  another 
elastic  modulus  other  than  Young's  modulus. 

The  free  torsional  vibration  of  a  cantilever  plate  composed  of  a  non- 
coated  substrate  is  considered.  The  objective  is  to  determine  the  shear 
modulus  of  the  film  from  the  shift  of  the  natural  frequency  from  that,  for  the 
non-coated  to  that  for  the  film-coated  substrate. 

Figure  25  shows  geometrical  symbols  for  the  cantilever  and  the  choice  of 
coordinates  with  the  x-axis  passing  through  the  centroid  of  each  cross 
section.  The  analysis  was  originally  carried  out  for  a  substrate  coated 
equally  on  the  top  and  bottom  surfaces  with  the  total  film  thickness  d^. 
However,  the  results  obtained  hold  equally  well  for  the  coated  substrate  with 
the  total  film  thickness,  regardless  of  the  thickness  ratio  between  the  films 
on  the  top  and  the  bottom  surface,  if  the  film-thickness  d^  is  much  smaller 
than  the  substrate-thickness  d. 
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Fig.  25. 


Geometry  and  Coordinates 
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•;  ir  j  I  vibration  of  a  beam  is  most  conveniently  described  by  the  ang^e 
of  me  cross  section.  For  a  cantilever  snowr.  on  Fig.  25 


d  -  b  <  <  i , 


a.  vi oration  is  governed  by  the  following  eigenvalue  problem: 
(GK  -  GfKf)  =  (oJ  +  ofJf)  ^ 

3  X  3 1 


9=0  at  x  =  0, 


—  -  0  at  x  :  !  (38; 

-»  Y 


s  without  subscript  refer  to  substrate,  those  with  subscript  f  tc 
shear  modulus,  o  the  density  and 


■  r  ~  ,  J  =  —  ( b  ♦  d  )  , 
0  ■  2 

bdp  2  -  bd 

K(,  -  (d  +  ddr  +  dfj  — 


od,  pp  i  go.. 

—  (b  ♦  3d  +  3dd.  +  dt.)  1  —  i: 

d  I  I  £ 


of  this  formulation  are  presented  in  the  Appendix  B. 

/,  let  p  and  p  denote  the  angular  frequency  for  the  fun 
mated  arid  the  uncoated  substrate,  respectively.  Further-, 


/  i_\  »\  /  U  X  \  ipt 

6 ( x , t )  -  A  sin  (—  -|  e 


ition  of  9  from  (40)  into  (37)  yields,  in  view  of  (39) 


n  -  1  |_ i _  G,  i/a 

po  '  4  1  ,2  ,2  o  * 

b  +  d 


(df  -  0) 


GK  +  G p  , 

—  I - Lii1' 

2S.  lPJ  +  DfJf' 


( df  *  0) 


Expand  the  right  hand  side  of  (40)  into  a  power  series  of  df/d  by  use  of 
;39)  and  (41)  and  keep  the  first  two  terms  to  get 


d„  Gr  ,2  _  .2  p„ 

r  ,  1  f  r  0  f  b  -v  3d  f  i 

Pol'+2  d  ^  G  “  2  2  p  > 

b  +  d 


Ap  =  p  -  pQ. 

Then,  the  shear  modulus  of  the  film  is  determined  from  (43)  in  terms  of  the 
frequency  shift  Af  as 

=  ‘w> 


which  is  the  desired  result. 
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A:u  END  1 X  A.  DERIVATION  OF  THERMAL  DISTRIBUTION  AND  DEFLECTION 


TEMPE KA  :'J RE  D I STR I  BUT  ION 


Temperature  distributions  in  the  substrate  during  and  after  the  coating 
deposition  are  well  represented  by  solutions  to  the  one-dimensional  heat  con¬ 
duction  problem  of  a  slab  bounded  by  planes,  z  -  0  and  d.  At  the  initial 
instant,  tne  temperature  T(z,t)  is  equal  to  the  ambient  temperature  I 
During  the  coating,  we  assume  that  a  constant  rate  of  heat  flux  qQ  per  unit 
area  is  supplied  at  the  front  surface  and  heat  radiates  at  the  rear  surface 
into  the  surrounding  medium  of  temperature  TQ  according  to  Newton's  law  with 
coefficient  of  surface  heat  transfer  H.  After  the  end  of  deposition,  it  is 
assumed  that  the  heat  radiates  from  the  front  and  rear  planes  with  a  linear 
initial  temperature  distribution. 

To  simplify  the  formulation  of  the  analysis,  we  introduce  the  following 
dimensionless  quantities: 

5  =  z/d 

•t  -  <t/d^ 

c  =  dH/K  (Al) 

vU,t)  =  [T(  z,  t )  -  T  Q  ]  /  [  T ,  -  T0] 

60  =  dvIHT,  -  T0)| 

where  K  and  <  stand  for  the  conductivity  and  diffusivity  of  the  substrate, 
respectively,  and  T1  denotes  the  equilibrium  temperature  at  the  front  surface 
when  the  deposition  is  continued  indefinitely. 

in  this  dimensionless  formulation,  the  heat  conduction  equation 


3  v  3  v 
2  =  3t 


(0  <  5  <  1, 


0  <  T  <  “1 


tht  deposition  problem, 


(0  V  ^  1), 


3t  t  :  0, 


L  V  —  0 

d 


at  [  =  1 ; 


and  a i n, i  lardy,  for  ’"he  post-deposition  problem, 


v  U,0)  =  a  + 
av 

— ^  -  £V  -  0 

H  P 

3  V 

— ^  +  ev  =  0 
3  C  P 


(0  '  j;  <  i), 


at  {  :  0, 


at  e  -  1 . 


The  terms  vd  and  vp  denote  the  solutions  to  the  deposition  and  the  post¬ 
deposition  problems,  respectively.  Expressions  representing  them  can  be 
obtained  by  such  standard  procedure^}  as  the  Laplace  transform: 


v3U,.l  =  e0;i  .  |  -  t)  -  HeB0  [ 


’ 

e  cos(X  £) 


n=1  X  [X  +  e(1  +  e)]  sin(2x  ) 
n  n  n 


where 


X  tan  X  =  e 
n  n 


-4<2t 

“  e  n  cos[  2c  ( (,  -  1/2)  ] 

Ma*§)  l— - 0 - 

n=1  [4k  4  e(2  +  e)]  cos  < 

n  n 


-4y2t 

-  e  n  sin [2y  U  -  1  1/2)] 

4(wf)b  l  - 5 - 2 - 

n  =  1  [ 4 y  +  e(2  +  e ) ]  sin  y 


wr.ere 


tan  k  -  — 
n  n  2 


Y  cot  v 
n  n 


Equations  (A5)  and  (A?)  involve  infinite  series.  Their  convergence  becomes 
slow  when  they  are  used  for  very  small  j  values.  Expressions  suitable  for 
small  t  values  are  established  by  asymptotic  inversion  of  the  transformed 
solution.  We  list  here  only  the  expression  for  v^ 


vdu’T)  =  2V 


ierfc  -  20o L(2  -  c,t) 


where 


-2BqL(2  +  ^ , t )  +  . . 


1  —  X 

ierfc(x)  =  ~TJ  2  e  ~  x  erf'c  x> 

TT 


L( t , t )  =  i  '  ierfc  ( — |  exp(e£  +  e2t )  erfc( — fyj  +  et'7')  ( A 1 0 ) 


-  -  erfc  f — f— 
£  2/' 


It  is  observed  from  Eq .  (A5)  that  the  steady  state  temperature  demands  a 

constrained  relation 


o  "  1  +  e  ’ 


(All) 


which  is  the  consequence  of  Eq .  (A1)  and  lim  vh(0,t)  -  1. 

“X  -*-ro  u 

In  ordinary  coating-depositions,  the  parameter  e  is  extremely  small.  For 
such  cases,  the  roots  and  Yn  of  transcendental  Eqs.  (A6)  and  (A8)  are 

closely  approximated  by 


fWY! iuf nArf WAkAjftLlT njftijf . ic'f il’TjOu aJCiJCbJl js  -  jA.  a  k  A 


X  =  ( n  -  1 )  n  +  t  [  (  n  —  1 )  t!  i  (  A 1 2 ) 

n 

(n  --  2,  3,  .  .  ■  ), 

-  (n  -  1)t  +  e/[2(n  -  1)n]  (  A 1 3 ) 

(n  =  2,  3, ■ • ■ ) , 

\  =  (2n  -  1)  ^  +  e/[2n  -  1)*]  (n=1,  2,...).  <  A 1 4 ) 


When  e  <  <  1,  simple  asymptotic  expressions  for  Vp  and  Vp  are  obtained  from 

Eqs.  (A5)  and  (A7),  respectively,  after  the  use  of  Eqs,  (All)  through  ( A 1 4 ) : 


vd(£;’T)  = 


1  -  e 


-ex 


-  eU  -  (1 +  f-) 


( A15) 


VpU.x)  =  (a  +  |)(1  ♦  y  h  -  c(t  -  |)?]  e'2eT 

Uh  .  2 

+  ~2  sin[  U  -  1/2)*]  e-11  T.  (  A 1 6 ) 


These  expressions  supply  satisfactory  values  except  for  very  small  x . 

II.  THERMAL  DEFLECTION 

Once  the  temperature  in  a  substrate  is  determined,  the  thermal  deflection 
can  be  predicted.  To  calculate  the  thermal  deflection,  it  is  assumed  that  a 
circular  wafer  of  radius  a,  thickness  d  and  coefficient  of  thermal  expan¬ 
sion  a  experiences  temperature  variation  T(z)  across  the  thickness  over  a  con¬ 
centric  circular  region  of  radius  b.  It  is  further  assumed  that  the  boundary 
is  simply  supported.  From  the  known  solution"^,  the  deflection  at  the  center 
is  obtained  as 


'  .  c.  u 

w(0)  -  [1  -  (1  +  v)  log(b/a )  ]  J  T(z)(z  -  d/2)  dz.  ( A 1 7 ) 

d3  0 

The  tip-deflection  of  a  cantilever  plate  of  length  l  with  the  temperature 
distribution  T(z)  across  the  thickness  is  known  to  be 

,  2  d 

-  -  —3-  J  T(z)(z  -  d/2)  dz . 


6 


{  A 1 8  ) 


In  view  of  Eq .  ( A 1 7 )  ,  the  integral  in  Eqs.  ( A 1 8 )  and  (18)  can  be  written  in 
terms  of  the  dimensionless  temperature  [ Eq .  ( A  1 )  ]  as 


d  P  1 

{  T(  z )  ( z  -  d/2)  dz  :  (T  -  T  )  d  J  V(c- OU  -  1/2)d*. 
0  0 


(  A 1 9 ) 


APPENDIX  B.  DERIVATION  OF  EQUATIONS  (37)  AND  (38) 


P 1 

We  follow  the  formulation  procedure  adopted  by  Rao  and  Gupta  for  the 
torsional  vibration  of  a  cantilever  plate  and  take  the  correction  due  to  the 
coating  film  into  account. 

Assume  the  displacement  field  as 

U  =  4>9  ,  U  =  -z0.  U  =  ye  ( B 1 ) 

where  <t>(y,z)  denotes  warping  of  the  cross  section  per  rate  of  twist  and  the 
subscript  x  preceded  by  a  comma  represents  differentiation  with  respect  to 
x.  This  rule  will  also  apply  to  all  space  and  time  variables.  Substitution 
of  displacements  from  Eq .  (B1)  into  the  strain-displacement  relation  leads  to 


xx 


,  xx 


e  *  c  =  0 

yy  zz 


e  =  ■=(  $  -  z)  e  , 

xy  2  , y  ,x’ 


E  :  0 

yz 


( B2 


e  z  -  {$  +  y ) 8 

xz  2  ,z  J  ,x 


Stresses  in  the  substrate  associated  with  strain  ( B2 )  are 
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xx 


E$e 


o  =o  =  o  -  0, 

yy  zz  yz 


0 

xy 


C(».y 


7  )  8 


0 

xz 


Stresses  in  the  film  will  be  obtained  from  (B3)  by  replacing  E  and  C  by  and 
Gj-;these  equations  will  be  referred  to  as  (B^3). 

The  boundary  conditions  for  the  cantilever  plate  are 


U 
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U 
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0, 
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xx 


0 


at  X  :  0 


a 

xy 


o 

xz 


0, 


o  =  0 

XX  ,  X 


at  x  :  l 


(BA) 


we  then  calculate  the  strain  energy  and  the  kinetic  energy  of  the  system 
and  appeal  to  Hamilton's  principle  to  arrive  at  the  governing  differential 
eeua t ion  : 


\  l  I  +  z.  r  i  ~ )  6  -  (.  oK  +  GrK  j  9 

ft  , XX XX  f  !  ,XX 


0T  +  ,  I  \  0 

2‘  f'  ,xxtt 


°J  +  Vr  e,tt  0 


02dA, 


K  :  r  r  [  t>  -  i )  +  ( $  +  y  h  '  d  A  , 

A*  ’ y  ’ z 


J  =  rr  (v" 


wh  :  o  the  expressions  for  1^,  K  _r,  to  be  called  as  ( B^7)  are  obtained  from 
IB?)  by  replacing  A  by  A,-  in  which 


y  ,  a  I  - 


y  *  A  ■ 


2  -  z  -  2 ; 


i  d 

:y.-7 !  - 


b  .  bi 

—  <.  7  <  — 

2  -  z  '  2' 


If  out  interest  is  to  find  the  fundamental  frequency,  we  may  assume  the 
mode  shape  as  given  by  Eq .  (42).  Substitution  of  6  from  Eq.  (40)  into 
Eq .  ( B6 i  shows  tnat  the  first  and  the  third  term  in  Eq.  (B6)  can  be  neglected 
"hen  b  •'  l.  Therefore,  the  original  problem  with  the  differential  Eq.  (B6) 
arid  the  boundary  conditions  Eq .  (B5)  can  be  replaced  by  a  simpler  problem  con¬ 
sisting  of  the  second  order  equation  (Eq.  37)  and  the  boundary  conditions  Eq. 
(  -;H ;  . 

For  thin  rectangular  cross  section  (d  ■  b),  the  warping  function  is 

approximately  given  by  (Ref.  22). 


Then,  expressions  for  K,  J,  and  are  evaluated  from  (B7),  (B^7)  and  (B8) 
as  g iven  in  ( 37) . 


The  analysis  in  the  Appendix  has  been  carried  out  for  the  case  where  the 
film  is  deposited  on  the  top  and  the  bottom  surface  of  the  substrate  by  an 
equal  thickness  d^/2.  On  the  other  hand,  the  final  results,  (41)  and  (42), 
refer  to  the  case  where  the  film  is  deposited  only  on  the  top  surface.  it  is 
to  be  pointed  out  here  that  these  results,  (41)  and  (42),  hold  true  regardless 
of  the  ratio  of  films  on  the  top  and  the  bottom  surfaces  since  the  effect  of 
such  ratio  appears  only  for  the  higher  order  terms  of  d^/d,  not  the  first  two 
terms  of  the  Taylor  series  expansion  in  (41). 
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